Diffractions of Bose-Einstein Condensate in Quantized Light Fields 
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We investigate the atomic diffractions of a Bose-Einstein condensate in quantized light fields. 
Situations in which the light fields are in number states or coherent states are studied theoreti- 
cally. Analytical derivation and numerical calculation are carried out to simulate the dynamics of 
the atomic motion. In condition that atoms are scattered by light in the number states with im- 
balanced photon number distribution, the atomic transitions between different momentum modes 
would sensitively depend on the transition order and the photon number distribution. The number- 
state-nature of the light fields modifies the period of atomic momentum oscillations and makes 
forward and backward atomic transitions unequal. For light fields in coherent states, no matter 
the intensities of the light fields are balanced or not, the atomic diffractions are symmetric and 
independent on the transition order. 

PACS numbers: 42.25.Fx, 42.50.Ct, 03.75.-b, 67.85.Hj 



I. INTRODUCTION 

The coherent interaction between matter and electro- 
magnetic fields within various kinds of physical system 
is providing a useful platform for developing concepts 
in quantum optics and atomic molecular physics. The 
experimental realization of Bose-Einstein condensation 
in dilute atomic gases [l|, greatly facilitates the investi- 
gations of the interactions between the atom and light 
on a macroscopic scale I n the past few years, 

atomic scattering process in optical fields has been inten- 
sively studied, such as matter- wave supperradiance 043 > 
quantum phase transition [l^, E| , cavity optomechanics 
[12l Il3|. raising a variety of striking discoveries. 

In physical systems involving light-atom interactions, 
various kinds of experimental configurations have been 
developed or proposed to achieve the desired atom-field 
coupling JEE 22 1 . Especially, by using high-quality res- 
onators 

PUT 

strong coupling regime of the light- atom 
can be reached, where atoms coherently exchange pho- 
tons with light fields. However, in the previous light- 
atom interacting models, the derivations usually use 
the classical treatment for the electromagnetic fields, 
in which the light fields are recognized as amplitude- 
modulated plane waves with slowly varying amplitudes; 
or even in a quantized treatment for light fields, mean- 
field approximations is often introduced for the light 
fields. Atomic diffractions in such situations have been 
well studied, however, in most conditions, the subtle ef- 
fects of the atomic motion induced by the quantum na- 
ture of light has always been ignored. 

In this paper, we study the diffractions of a Bose- 
Einstein condensate (BEC) in multimode optical fields, 
with a full quantum treatment of the light fields. In the 
cases that the quantized light fields are in different kinds 
of states, exotic behaviors are induced in the atomic mo- 
tion by the quantum-nature of light. This paper is orga- 



nized as following, in section II, we introduce the physical 
system and develop the theoretical model of the atomic 
diffractions in a full quantum treatment of the light fields. 
In section III, we analyze how the quantum status of the 
pump fields can modify the atomic motion in the cases 
that the optical fields are in number states. In section 
IV, we turn to the case where the light fields are in co- 
herent states. In section V, we arrive at a summary of 
our results. 



II. THE THEORETICAL MODEL FOR ATOMIC 
DIFFRACTION IN QUANTIZED LIGHT FIELDS 

The physical system under study is illustrated in Fig.[TJ 
where an elongated BEC is trapped along the horizon- 
tal axis of a triangle ring cavity. The cavity consists 
of three mirrors which are placed in a way such that 
the light reflected inside the cavity forms a closed loop. 
There are two mutually counterpropagating modes, i.e., 
the clockwise and anti-clockwise modes, coupled equally 
to the atoms. The light fields in these two modes 
are also referred as the forward (or right-going, with 
wave-vector/cfl) and backward (or left-going, with wave- 
vector kr,) ones according to the geometrical configura- 
tion. These two cavity modes are degenerated in fre- 
quency, ujji — cut, — u) c , so the wave- vectors of the two 
cavity mode satisfies k^ = —kj,. The atoms within the 
condensate are recognized as two-level atoms with one 
internal ground state \gi n ) and excited state |e, n ) repre- 
senting the internal degrees of freedom for each atom. In 
the large detuning limit, the upper internal state \ei n ) 
of the atoms can be adiabatically eliminated [§, [13, Ell ■ 
For illustration, the condensate is chosen to be prepared 
with 87 Rb atoms, and the wave length of the light fields 
is A = 780nm with a detuning A = lo c — io a = — 1.5GHz 
of the £>2 transition M, [U [53 ■ The Hamiltonian of the 
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FIG. 1: (Color online) A schematic diagram of a two-mode tri- 
angle ring cavity containing an elongated Bose-Einstein con- 
densate trapped along one of its cavity axis. The damping 
rate k, of the cavity modes is small compared to the atomic 
motion inside the cavity. There are two degenerate light 
modes propagating in the clockwise and anticlockwise direc- 
tions inside the cavity. 



light-atom coupling system can be written as, 

r * r ft 2 ~ i ~ 

,ye= I dW(r) y - ^-V 2 + E(r) • D(r)J *(r) 

+ iy , d 3 rd 3 r*t(r)*t( r y mter ( r „ r ')*(r')*(r), (1) 

where ^(r), ^(r) are the creation and annihilation field 
operators for the atomic field, respectively. E(r) is the 
electric field operator of the electromagnetic field, and 
D(r) is the atomic dipole moment operator. Vi n t er (r — 
r ) represents the interatomic potential. The counter- 
propagating light forms a standing wave in the triangle 
ring cavity. The atoms in the standing light feel a peri- 
odical potential. The condensate is transversely tightly 
bounded, therefore its transverse freedom is frozen, and 
it is reasonable to recognize the condensate as quasi- 
one-dimensional ensemble in its longitudinal direction 
19]. This suggests that the matter field operators can 
be expanded with a set of Bloch functions, which are 
amplitude-modulated plane waves, 



*(x) = ^fe„e jk " x u(x), 



e- jk " x u*(x), 



(2) 



where u(x) is a periodical function along the direction 
the horizontal axis of the cavity, u(x + 4) = u(x), 
with the period one-half of the wave-length of the light 
wave. 6^ and b n are the creation and annihilation op- 
erators of the bosonic atoms in the nth Bloch mode. 
k n = 2nkfl is the wave vector of the modulated plane 
wave. For the absorption and subsequently stimulated 
emission of photons from the two modes of optical fields, 
the atoms would acquire a net momentum of even num- 
ber multiples of photon momentum, InKk^. Assuming 



the atoms are initially in the stationary state, after the 
scattering process, they will be excited to the higher 
adjacent momentum modes, which can be described as 
( e 2»k R x _|_ e - 2lk R x ) Uo (x), with Uo(x) the envelop of the 
stationary wave packet. If the atoms get further scat- 
tered, some of them may go to the even higher mo- 
mentum mode. For simplicity, we can approximate the 
amplitude-modulated plane waves with the wave func- 
tions of the free particles, and the wave functions of the 
atomic field now read, 



n 

=$>;<p„|. (3) 



In our case, due to the diluteness of the atomic gas, we 
are going to neglect the modification of wave functions 
from the atomic interation. Here, V is the volume of 
the condensate, (x\p n ) = ^j=exp(Jrp n • x) is the wave 
function of free particles with momentum p n — 2nhkR 
in coordinate space. The particle number of the atoms is 
normalized to one. At the same time, we have removed 
the hat from each operator of the matter field, using a 
mean- field treatment for the atomic field. 

We approximate that the discrete modes of the atomic 
field makes a complete set describing dynamics of the 
condensate. Plugging Eq. ^ into Eq. ([T]), we arrive at 
the following expression for the Hamiltonian of the light- 
atom coupling system, 

oo 

fi^m \Prn) (Pm| + hG 

in— — oo 
oo 

Y\ a k L ak B LP") (Pn-l\ + a l R a k L \Pn) (Pn+l|) , (4) 



n— — oo 



where we have adopted the full quantized version for 
the electric component of the light fields, E(r) = 

X^k<£k(ak + &k)>k = ^L,R, m which £k is the unit 



polarization vector; = ^Jhio c /2eV c has the dimen- 
sion of an electric field with V c the quantization vol- 
ume of the macro-cavity [2(| . ak and a k are the cre- 
ation and annihilation operators of photons in the corre- 
sponding light mode, respectively. fvuj m — 2(ftmfcz,) 2 /M 
is the energy of the mth unperturbed discretized mode, 
and M is the mass of a single atom within the conden- 
sate. G = — fl 2 / A is the effective two-photon atom-field 
coupling constant, with fl = — d e ^' ?k the single-photon 
atom-field coupling constant and d the electric dipole 
moment of the atom. The form of G comes from the 
adiabatic elimination of the internal excited state of the 
two-level atom, and this is legitimated by large detuning 
of the light fields. This Hamiltonian is rather reminis- 
cent of an ordinary lattice model in solid state physics 
in spite of an atomic interacting term. The second hop- 
ping term disturbs the motion of the free atoms and in- 
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duces exchanging of particles between different diffrac- 
tion orders. Such atomic diffractions are categorized 
into Raman-Nath diffractions with short pulse fields and 
Bragg diffractions with longer time of light-atom inter- 
action. 

To investigate the atomic motion of the system, it is 
usually convenient to work in the interaction picture by 
seperating the Hamiltonian (jl]) into the free part and the 
atom-field interacting part, then one has, 

oo 

f(t) = hG {4 L ^n\Pn) ( Pn -l\e- iS -W 

71 — — OO 

+ a{ R a kL \p n )(p n+1 \e- U +^}, (5) 

where 6± = ±l6iris r (n ± i), and v T = hk\/(AnM) is 
the atomic one-photon recoil frequency (v r ~ 3.77 kHz 
according to the parameters chosen in this system). 

The optical freedom of this system can be integrated 
out to extract the effective information of atomic dynam- 
ics. Unlike the classical treatment of the optical fields in 
which the optical fields are recognized as plane waves 
with slowly varying amplitudes, the atomic motion of 
the system can be explicitly modified by the status of 
the quantized light fields, not just by the strength of the 
pump beams [III, [Hf. To see this, we shall assume, for 
example, that the light fields are in number states. The 
Hamiltonian ([5]) conserves the total photon number of 
the system. Thus the wave function of the matter-wave 
condensate can be written as a linear combination of dif- 
ferent momentum states with the corresponding photon 
distributions among the two light fields, 

|* (<)) = Y *« \i> (n, N R ~n,N L + n, t)) , 

n 

where n = 0,±1,±2, refers to the diffraction order. Nr 
and Nl are the initial photon number distribution among 
the two light modes. \ip (n, Nr — n, Nl + n, t)) is the 
product of the wave function of the atomic center-of-mass 
motion and that of the optical fields, \p n ) ® \ipf)- In this 
stage, \ipf) is choosen to be the number state noted as 
\n\,n2) corresponding to that there are ri\ photons in 
the right-going light mode and photons in the left- 
going mode. Since the condensate is initially prepared in 
the stationary state, each atom needs to absorb n pho- 
tons from the right-going light mode in order to hop to 
the nth momentum mode. Therefore, the photon num- 
ber distribution of the two light modes is directly related 
to the atomic diffraction order. The equation of atomic 
motion now reads, 

ih | |* (t)) = r(t) |*(t)) 

oo 

= hG J2 {<i L a kR \ Pn ) (p n -i\e- i5 -V 

71 — — OO 

+a\ R a kL \p n )(p n+1 \e- iS +W}\*(t)). 



Keeping the photon degree of freedom, the above equa- 
tion can be further simplified to, 

*n(*) \n R ~n,n L +n) 
= a^a fcn vl/ jl _ ie - 4 ' 5 -* \ Ur - n + l,n L + n - 1) 

+al R a kL ^ n+1 e- l5+t \n R - n - 1, n L + n + 1) , 

Until now we have kept quantum character of the optical 
fields. Integrating the photon degrees of freedom by ex- 
plicitly executing the operations of the photon creation 
and annihilation operators on the corresponding number 
states of the light fields, we arrive at the following equa- 
tion of the atomic motion, 

*»(*) = ■ *„-i e - 15 - 4 + W£ ■ * n+1 e- i5 + t , (6) 

where W£ R are the forward and backward transition 
weights defined as, W R — —iG\J (ur — n + 1) {ul + n) 
and W£ — —iGyJ (hr — n) (nj, + n + 1). This equation 
has included the quantum effect of optical fields upon 
the atomic motion. The complicated expressions for the 
transition weights are resulted from the property of num- 
ber state of the optical fields that when it is operated by 
the corresponding creation and annihilation operators, 
we have (r\n) = \Jn + l\n + 1) and a\n) — y/n\n — 1). It 
can be directly observed from Eq. © that, the transition 
weights depend not only on the diffraction order but also 
on the photon distribution among the two light modes. 
Thus, unlike the classical case of plane wave approxi- 
mation for the light fields, there is no definite equality 
between W% and W£; \W£\ > \W£\ and \W£\ < \W£\ 
are both possible, depending on whether Nr < Nl or 
Nl > Nr. Table (HJ) lists out the transition weights of the 
first few order atomic scattering processes for Nr = 80 
and Nl = 10. It shows that for the first few order scat- 
tering processes, \W R \ < \W£\ are always satisfied. 



Scattering Process 


\WR\ 


\W£\ 


\WR ■ W£\ 


bo) -> \pi), \p-i) 


28.46 


29.66 


844.12 


\Pl) |Po), |P2> 


29.66 


30.78 


912.93 


b-i) -> bo), b-2) 


27.16 


28.46 


772.97 



TABLE I: The transition weights in dependence on the 
diffraction order in unit of atom-field coupling constant G, for 
the initial photon number distribution, Nr — 80, Nl = 10. 



III. LIGHT FIELDS IN NUMBER STATES 
A. Raman-Nath regime 

To solve Eq. © in the full time-domain, one has to 
appeal to numerical calculations. Fortunately, in the 
Raman-Nath regime where the atoms interact with the 
light fields for rather short time, the low order diffrac- 
tions dominant the whole atomic scattering process, and 
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FIG. 2: (Color online). Atomic probability distribution of 
the first few discrete momentum modes in the limit of short 
light-atom interacting time, for an initial photon number dis- 
tribution, Nr — 80, Nl = 10. Black solid curve refers to 
n = mode; red solid curve , n = 1; blue dash-dotted curve, 
n = — 1; green-solid curve, n = 2; wine dashed curve n = — 2; 
pink solid curve, n = 3; orange solid curve, n = —3. The 
light-atom coupling constant is chosen G = 0.7v r , with v T the 
atomic one-photon recoil frequency. 

Eq. ijJJJl can be solved analytically. In this limit t 0, 
the probability equation of the atomic motion takes the 
following form, 

*n(t) = ^fl ' *n-l + W£ ■ *n+l- (7) 

For macroscopic occupations of photons in the two light 
modes, the diffraction order n is small compared to the 
initial photon number Nr and Nl, i.e., Nr.l 3> h, we 
have Wjj ~ W£ ~ s/NrNl which can be observed in 
Table ©. Thus Eq. © has the following solution, 

W n 

^n(t) = i n e i9 (^J n m (8) 



where £ = 2 y / (W R l W£), J n (($) is the nth order Bessel 
function; 8 is an arbitrary phase angle. The popula- 
tion of atoms occupying the nth momentum state \p n ) is 
Pn {t) = |^„(t)| 2 = (^l) B [J n (#)] 2 . For the inequality of 
Wo and as show in Table (JTJ) , the atomic diffraction 
will display exotic behaviors. Fig.f2]shows the probability 
amplitudes of the atoms in different momentum modes 
evolving with the light- atom interacting time. It can be 
seen that, besides that the atomic populations in differ- 
ent momentum modes oscillate with time, the curves of 
the momentum oscillations of the ±nth modes do not 
overlap. Unlike the situations illustrated in [HHHI, the 



probability amplitudes of these ±nth momentum modes 
do not reach the minimum or maximum values exactly at 
the same time. For example in Fig. [2 there is a shift be- 
tween the first minima of the probability amplitudes for 
the ±lst order atomic momentum modes. Even though 
to study the atomic motion for long-period light-atom 
interaction from the data derived in this regime may 
not seem adequate enough, it implies that the atomic 
diffractions would be asymmetric about the zeroth order, 
and such asymmetries are induced by the number-state- 
nature of the optical fields. More detailed discussion and 
the accuracy of the analytical solution in this short time 
regime will be demonstrated in the later part of this pa- 
per. 



B. Bragg regime 

In the Bragg diffraction regime, where the light and 
atoms interact with each other for a longer time, the ex- 
ponential factors in Eq. ^ greatly modify the atomic 
motion, compared to the situations of the short time 
light-atom interactions (i.e., Raman-Nath regime) . How- 
ever, in this regime, the equation of atomic motion can- 
not be explicitly solved analytically. To obtain the in- 
formation of the atomic motion in the full-time domain, 
numerical calculation is carried out with a cutoff at the 
the n = ±10th momentum modes in order to maintain 
the accuracy. In fact, the validity of such a cutoff has 
been implied in the case of short-time light-atom inter- 
actions; for the very high diffraction orders, the atomic 
occupations in the relevant modes are always vanishingly 
small. 

Because of the dependence of the transition weights 
upon the diffraction order and the photon number distri- 
bution, numerical simulation displays interesting behav- 
iors of the atomic motion in the quantized light fields. It 
seems that the atoms can see the imbalanced intensities 
of optical fields, and thus take actions accordingly. Fig. [3] 
shows the probability amplitudes of atoms in the first few 
order momentum modes for different photon number dis- 
tributions. It is found that, the quasi-period of momen- 
tum oscillations is changing with respect to the diffrac- 
tion order. Fig. EJa) shows the probability amplitudes 
of the scattered atoms at the photon number distribu- 
tion Nr = 80, Nl = 10. Even for atomic diffractions of 
the same order, see the n = ±1 curves, the atomic pop- 
ulations oscillate at different frequencies. At t « 28/xs, 
the curves for the n = ±1 modes get to their lowest 
points, but if one goes to the details in this interval, 
it is found that there is a shift between the minima of 
two curves. That is to say the populations of atoms in 
these two modes are in fact oscillating with different pe- 
riods. Additionally, it is found that the atoms would fa- 
vor the forward and backward diffractions alternatively; 
the n = 1 order is sometimes stronger and sometimes 
weaker than the n = — 1 order, implying that asymmetry 
appears in the matter-wave diffractions. This peculiar 
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FIG. 3: (Color online). Atomic probability distribution on 
the first few discrete momentum modes in the regime of long 
light-atom interacting time, for the initial photon number dis- 
tribution, (a) N R = 80, Nl = 10; (b) N R = 80, N L = 10. 
[Line styles are the same with that in Fig [2^ a)] 

phenomenon is resulted from the quantum nature of the 
optical fields. The asymmetric atomic diffractions can be 
manipulated by adjusting the imbalance of the photon 
number distributions in the two light modes. Fig. |3jb) 
shows the atomic distribution of the diffracted atoms for 
Nit = 200, Nl = 20. Similar diffraction behavior is ob- 
served as that shown in Fig. Ela). However, the atomic 
distributions in this case are oscillating at higher frequen- 
cies. Thus the imbalanced photon number distribution 
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FIG. 4: (Color online), (a) Atomic probability distribution 
of the first few discrete momentum modes from numerical 
calculation in the Raman-Nath regime, for the photon number 
distribution Nr = 80, Nl = 10. [Line styles are the same with 
that in Fig Ufa)], (b) The error of the analytical solution 
obtained in the short time limit of the light-atom interaction, 
compared to that obtained in the Bragg regime. 



of the light fields modifies the dynamics of the atomic 
diffraction, and the order-dependence of the transition 
weights induced by the number-state-nature of the light 
fields makes the atomic momentum oscillation very irreg- 
ular. 

The above result is quite different from that obtained 
from the traditional analysis, where mean-field approxi- 
mations are applied to both the atomic field and optical 
fields. The photon creation and annihilation operators 
are directly replaced with their mean-field values. Even 
though this is efficient for extracting qualitative infor- 
mation about the atomic dynamics, it ignores the details 
of influences from the quantum nature of the light fields 
upon the atomic motion. 

To check the accuracy of solution derived in the short 
time limit of the light-atom interactions in the last sec- 
tion, numerical simulation for Eq. is carried out for 
the photon number distribution Nr = 80, Nl = 10, and 
the result is shown in Fig. EJa). It is found that the 
data obtained here is consistent with that obtained in 
the Bragg regime. There are displacements between the 
minima (or maxima) of the conjugate ±n curves, and 
the probability amplitude of the n — 1 momentum mode 
can be either larger or smaller than that of the n = — 1 
mode. This is similar to the result displayed in Fig.[3Ja). 
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However, we have to emphasize that, using Eq. (J7J to 
describe atomic motion is only valid in the Raman-Nath 
regime with short time light-atom interactions; thus data 
in Fig. IDJa) for large t is not trustful. Fig. 0Jb) shows 
the error of the data in this regime compared to that 
from Fig. EJa). The accuracy function D a is defined as 

£>"(*) = |fr||2 + |^||2 , where is the probability am- 
plitude of the atomic motion in the Raman-Nath regime 
obtained in section II. A, and is that obtained in the 
Bragg regime. It is found that, for short time, t < 5/xs, 
the error of the analytical solution Eq. © is less than 
1%. So it is a good approximation of the atomic motion 
in this regime. 

IV. LIGHT FIELDS IN COHERENT STATES 

Besides that the properties of atomic diffractions are 
sensitively dependent on the intensities of the pump 
fields, the status of the light fields can also impose im- 
plicit corrections to the atomic diffractions. To make 
comparison to the case that the lights are in number 
states, we turn to the case that the two optical fields 
are in coherent sates \ocr) and \ocl), respectively, an 
and chl are the eigenvalues of the annihilation operators 
of the two light modes on the relevant coherent states. 
The equation of the atomic motion reduces to 

* n (t) - W co ■ * n _ ie - tf -* + W* co ■ * n+1 e- i5 + t , (9) 

where W co and W* a are conjugate with each other, W co = 
ana* L . Generally, or and ai, are complex numbers. If 
an and a_L are chosen to be real, then W co = W* a = 
a roil and the eq uation is of the same form as that ob- 
tained in [2H l22j |. By numerical simulation, the infor- 
mation of the atomic motion is obtained in the full-time 
domain. Fig. [5] shows the probability amplitudes of the 
first few momentum modes at W co — 60 + 40j. It is found 
that the curves of the ±nth orders coincide with each 
other; thus the forward and backward atomic transitions 
are equally favored, and no asymmetry appears in the 
atomic diffraction. This can be seen from the fact that 
the absolute values of transition weights for the forward 
and backward transitions are the same, and their phase 
difference can be absorbed in the definition of the phase 
factor 5± as constant, thus the equation of the atomic 
motion still takes the symmetric form as that in [2ll . l22j . 

In this case, the transition weights are independent 
upon the diffraction order. The data obtained here 
resembles that with classical treatments of the optical 
fields, where the electromagnetic field are recognized as 
plane waves. Comparing the results with that in the 



case where the optical fields are in number states, even 
though the equation of the atomic motion is also derived 
in a full quantum treatment of the light fields, only sym- 
metric atomic diffraction can be obtained here; this is 
determined by the specific status of quantized light fields. 
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FIG. 5: (Color online). Atomic probability distribution of 
the first few discrete momentum modes in the condition that 
the light fields are in coherent sates, \an) and \ai,). The 
parameters are chosen such that W co = Q-rQ-l — 60 + 40i. 
[Line styles are the same with that in Fig[2ja)]. 



V. SUMMARY 

In this paper, we investigate the atomic diffractions 
in quantized light fields. The quantum nature of the 
light fields induces exotic phenomena in the atomic scat- 
tering processes. For light fields in number states, the 
atomic transitions among different discrete momentum 
modes depend on both the transition order and the ini- 
tial photon number distribution. The quasi-periods of 
the atomic momentum oscillations are implicitly modi- 
fied, and the atomic diffraction is no longer symmetric. 
However, for light fields in coherent states, the atomic 
diffractions are still symmetric, this verifies the coherent- 
state-approximation of the light fields which is frequently 
utilized in the cavity quantum electrodynamics. Situa- 
tions of other kind of states of quantized light fields, such 
as arbitrary linear combination of number states, can also 
be investigated using the methods introduced in this pa- 
per. The results obtained in this paper can be extended 
to other light- atom interacting systems for enabling new 
physics in atom optics. 
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